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ABSTRACT 


A study of a multimode third order instrument servomechanism is 
presented, Particular emphasis is placed on the effects of the initial 
condition of a compensating network, switched in as the system enters 
a linear operating mode, upon the behaviour of the output. The optimum 
values for the initial condition as a function of velocity and accelera- 
tion of the servomotor is derived, using the integral squared error as 


a criteria. The results are confirmed with a computer simulation. 
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I, INTRODUCTION 


Multimode control systems are those which are operated under two 
or more controlling laws which are in turn switched on or off according 
to some pre-established rule. Here, a particular case is studied, in 
which an instrument-servo is used to position a device in two main 
steps: 1.- First the motor is driven open loop at full forward speed, 
and then slowed down by reversing its voltage — both without feedback; 
2.- Finally, within a given error limit, a feedback loop is closed, 
zeroing the error according to the linear-system response, 

It is the second of these steps that is studied with greater detail 
as it is here where the optimization technique of this work applies, 

In general, in this type of system a third order type one instrument 
servo 1s used which is unstable under the desired gain conditions. 
Therefore, a compensating network has to be added to the forward path 
of the loop which is a lead network in this case. In practice, the 
initial condition on the capacitor of this lead network is usually made 
equal to zero, as it is easier to have a short circuit across the capac- 
itor than to provide an extra voltage source for this initial condition. 
No investigation has been found in the literature as to the effect of 
this initial condition on the time response of this type of control, and 
it is the purpose of this study to determine its effect on the servo 
performance. 

To accomplish this, a typical third order instrument servo of the 
linear actuator type has been modeled and the whole problem simulated 
with the aid of the system/ 360 Continuous System Modeling Program of 
IBM ( S/360 CSMP ) in the IBM 360/67 Computer. 








o> 


The results obtained from this special case simulation will be 
discussed in the light of producing a more general criteria for similar 
cases. 

The performance criteria to be used will be that of the Integral 
Squared Error, since this function provides the necessary information 
on minimum error which is desired here. Minimum settling time will be 


observed from the data obtained from the computer simulation. 


II. BASIC CONCEPTS ON DISCONTINUGUSEsy Sms 


Consider first a general block diagram of a discontinuous control 
system as shown on Fig. 2.1. The general control is achieved by the 
control computer which can be an analog or digital computer, an ampli- 
fier with relay, or, a simple relay circuit. The controlling signals 
may come from any point in the system, the input and output being the 
most frequently used. The control computer uses its input signals to 
provide the desired outputs which control the parameters of the plant, 
the controller and amplifier, or the feedback signals. 

From the analysis of the input signals, the input signal analyzer 
provides the required signals for the control computer. This analyzer 
may consist of a differentiator or integrator network, a sampler, a 
potentiometer or a function generator. 

The controller and amplifier network may consist of an amplifier 
with relay with R-C network, or a low power amplifier circuit which 
provides a compensated signal for the plant drive. The control computer 
can control the gain of the amplifier and the compensating circuit. 

The plant is the object to be controlled by the rest of the devices. 
In electrical control systems it may consist of a motor driving a load, 
its shaft being designed to follow the input signal. The plant must also 
be designed to meet the specifications of the system, which in the given 
case would be settling time, overshoot, transient response, etc. The 
control computer can vary the gain, the damping and/or the inertia of 
this plant. 

The feedback block consists of one or many signals which are a 


function of the outputs. The design will consider the number of signals 
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to be used, and the control computer can control the polarity, instant 
and period for them to be applied, as well as their gains and character. 

From this general structure, many special cases can be derived by 
placing limitations or specifications on the blocks indicated above. 

By using the controller and amplifier block and the plant, various 
control methods have been developed, The basic theory along these lines 
is to design the controller to make the switch operate along a curve in 
the phase plane, or a surface in the phase space that is called the 
optimum switching curve or hypersurface (Ref eae 

For the continuous operating mode of the system, a compensating 
network is used, and if blocks 2 and 3 are used, the compensation is the 
cascaded type. This form of compensation seems predominant in this 
field, although other types of compensation could be used. The theory 
of cascaded lead - lag networks has been largely discussed and developed, 
the Bode diagram being used with great advantages in the design. 

The reason for using discontinuous systems arises from the need of 
applying a limited amount of power to accomplish traveling in minimum 
time. These types of systems have been proven to be the best form of 
control even under economic considerations. 

Discontinuous systems are based on the theory that it is best to 
use maximum drive at the beginning of the trajectory, and then brake 
the system at a proper point. Thus, full voltage is applied at the 
beginning, and full reverse at the switching point. With the disconti- 
nuity, the character of the system is entirely different when the state 
point is at different sides of the hypersurface. The decaying part Oa 
the trajectory should be made as close to the E vs E plane as possible 


Rete 24. 
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Comparing the continuous system with the discontinuous, the phase 
plane plot of E vs E serves best to give a better insight of the differ- 
ence, Fig. 2.2 shows a typical phase plane of a discontinuous system, 
whilst Fig. 2.3 presents that of a second order continuous system. 

The continuous system is normally designed to have some overshoot 
in the trajectory, because an operation in the overdamped or even 
critical condition provides a slow response. In the discontinuous 
system the acceleration trajectory can be made unstable, and the die 
out trajactory a straight line near the E axis in the phase plane 
(fast eigenvector). 

This general discussion has only concerned itself with possible ways 
to use discontinuous theory to improve the response character of systems. 
The general switching theory and response character will not be discussed 
here in more detail. The main concern in later chapters will be the 
effect of entering a final phase which will be continuous; the use of 
a lead type conpensating network, and the influence of the initial 


conditions, in particular that of the compensating network. 


1] 





Mme 


Fig. 2,2 Typical E vs E Plot for a Discontinuous Operating System 








Fig, 2.3 Typical E vs E Plot for a Second Order System 





III. PLANNING THE EXPERIMENT 


A. ANALYTIC DISCUSSION 

In order to study a case of how the initial conditions on the lead 
network would affect the performance of a system, the problem was simu- 
lated using the characteristics of a linear actuator, such as the ones 
used in disc memories for computers to move the arm which carries the 
magnetic head. 

The system is supposed to move different eee the only Known 
variables being the initial error and the actual position of the head 
at all times. Minimum time and final error are of the essence, thus 
leading to a discontinuous bang-bang type of operation. The final part 
of the trajectory, i.e.: when the head reaches the edge of the track 
to which it was ordered, is accomplished in a closed loop mode. An 
extra gain is introduced to the system, and thus a filter has to be 
added in cascade in order to achieve stability, the filter being a lead 
network. 

Once a switching is determined, in order to have a real starting 
condition of velocity and acceleration of the motor as it enters the 
closed loop mode, a computer simulation using the S-CSMP/360 was 
implemented to study this final phase. 

1. Equations of the Physical System 

The known characteristics of the motor in practical MKS units 


Were. 


Newton 


Force constant = Kk. wee Ampere 


Motor inductance = L, = .01 [Henry ] 





Motor resistance = R, = 5 [Ohms ] 


Mass of motor = M = 0.45 [Kg.] 


ees Newton:sec 

Viscous damping = f = 17.9 Meter 

: Volts +sec 
15.6 Meter 


Back EMF constant = K 


b 


Armature voltage = v, [Volts ] 


Mx + TX = Nala (3.918) 


Kx + ale Rede ae (cee 


which, using Laplace transforms and substituting for i,, gives 


( 2 Re Sta 2 
sM+ sf) ot t+ ~~ (s Mt sf) + sk, X(s) 





V(s) = 2" c 
m 
or 
G(s) = Oy ee eS oe 
VS) gp te gg (ee eet 
S K a i b 


which is the transfer function for the motor with voltage source feed. 
Similarly, the study could be made by considering a current 


source feed, in which case the transfer function would be 


>< 


Gy n (3.4) 


ES) Ssh ct ) 








Going back to the voltage feed transfer function, and substi- 


tuting the given values of the motor characteristics, the open loop 





gain is: 


(Gi (3.5) 


s[s* + 538.65 + 72741] 
From this expression, and by simple inspection, it is known 


that 


w, = 269.7 a 


0) epeisis 


JN 
Il 


for the quadratic factor in the denominator. 

Thus the servomechanism is a type one, third order system, with 
two complex poles almost critically damped. The Bode diagram for such 
a motor is shown in Fig. 3.1. 

From the Bode diagram, it is seen that the KY of the motor is of 
the order of 0.047, which will represent an exceedingly slow system for 
the final closed loop phase of the operation, as compared with the desired 
settling time. 

2. Filter Design 

By further inspection of the Bode diagram, it is determined that 
an amplifier which should provide an extra gain of 55,000 ought to be 
used. This new characteristic function gives an unstable closed loop 
system, which in turn leads to the necessity of adding the lead network. 


A general and most simple form of lead network is represented by 


g(s) = a(st + 1) (ena) 


ast + | 


and its physical implementation is shown in Fig. 3.2, where % is the 


relationship between the corner frequencies of the filter as referred 
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to the Bode plot, i.e.: 





In order to obtain the maximum phase correction with the filter, 
an a of 0,1 is used [3], and by trial and error process the system 
is compensated to provide the Bode plot shown in Fig. 3.4, the filter 
having the following characteristics 
a= ,] 
t = 0.004 sec 
Thus 


G(s) = 0.1(.004s + 1) (a7) 


.004s + | 
The compensated system has now a transfer function 


a: G(s) = 5500(0.004s + 1) 3432.3 


a ae ek eT (3.8) 
0.0004s + 1] s(s_ + 538.65 + 72741) 


It is a stable system with a phase margin of 52 degrees, its 


foot iocus plot has been presented in Fig. 3.5. 


B. OBSERVATION OF THE EFFECT OF THE INITIAL CONDITIONS 

With the help of the flexibility of CSMP, several values of initial 
conditions for the filter can be implemented in one computer run, thus 
allowing a gross estimation of the direction of the variation of the 
performance with the said initial condition. 

The performance criterion used is that of the integral squared error 
in the linear mode, The optimum value of Vc is calculated to minimize 


this index. Since there are no weighting factors to consider in this 
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Root-Locus Plot for the 


xscale = 500 units/inch 


yscale = 5 units/inch 
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Compensated System 
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case, the problem is a very simple one, and the general expression for 


the optimal control 


t 
Jo = WUx(te) ste) + s¥ g(x(t),u(t).t) dt (3.9) 
0 


where to and te are the initial and final time; h and g are scalar 
functions; x(t) is the trajectory and u(t) the control. [Ref. 4] 


reduces to simply 


tf 
J= he bode (3.10) 
0 
in which 
Be gion) ule) eit) (32110 


and 1s given by the transfer function and solved by the program itself. 
Although flexible, CSMP places a number of restrictions on the use 
of subprograms and subroutines, therefore making it quite difficult to 
use the normal and known procedures for minimization and optimization. 
This leads to a rather crude method of determining the best suited 
initial condition for the filter, namely that of analyzing an ample 
range of values which show how the response of the system, and in 
particular the integral squared error, vary with the variable in 
question. Once a relative and feasible range of values for the initial 
condition is found, the best fit can be determined without much effort. 
1, Algebraic Development 
The solution for x(t) will now be pursued using the state 
variable form approach. 
From the motor equations (3,1) and 3.2), it follows that: 


X(t) = -ax(t) = bx(t) + ev, (t) (3.12) 


Za 











where a= ap a: a = 72741) 


Also, from the filter circuit in Fig. och 





i(t) = ip(t) + ip(t) (3.13) 
ie) Ge) 
et tee eee tt) (3.14) 
Ry R 
e,(t) = e,(t) + vo (t) (35 15)) 
Thus 
> lion i ail) 
ve(t) = =( RC” RC ) ve(t) + — (3.16) 
Rie 
and since RC = 7 
and R = QR, 
vo(t) = os le Vo (t) + Ze. (t) (sa) 


This result holds under the assumption that the motor does not 
load the filter circuit, which can be achieved in practice with the 
use of an isolating amplifier. With this in mind, the analysis can be 
furthermore pursued to solve for the relationship between the output 


x(t) and the initial condition on the capacitor V_ (0). 


Since the system is closed loop with zero input, the general 


expression 
E=R-C 
reduces to 
E=-C 
and therefore 
e,(t) = Ke(t) = - K x(t) (Sie) 


where K is the gain of the forward feed amplifier. 


From (3.17) and (3.18) 


v(t) = - PO ve(t) A 28 x(t) (3.19) 


Similarly, taking equations (3.12) and 3.15), and considering 
that 


e, (t) = v(t) 
then: 


x(t) = <ckx(t) © ax(t) - bx(t) <- ev, (t) (3,20) 


Expressing equations 3.19 and 3.20 in their state 


variable form, and letting 


x = yale) 
me Ce) 
Xq = x(t) 
Xq = Ve (t) 
then 
X(t) = Ax(t) + Bu(t) (3.21) 


24 


where 


0 ] 0 0 
0 0 ] 0 
AS and u(t) = 0 
=CK -a -b -C 
9K 0 0  -10 
Ae Af 


Using Laplace transforms, the solution to equation (3.21) is 


eS 


x(t) = £ [st - a}7!x(o) 


from which, the desired solution is the expression for X= Xe) 


enact 1S 


i ae b 0 
x, ( ) | 5(52 = re ‘i 10) ve : (s + ~)(s ct; Sta) x, ( ) 


t 


+ (stb)(s + *9)x5(0) + (s + *)x,(0) + (0 


mt) =fo' ee ae: (s + 10) (52 bsta)x(0) 
(stm)(stn)[(stp)* + 97] ; 


+ (stb)(s + 2)x(0) + (s + *2)x(0) + e() 


where 248.019 


Ss 


CSS 206 


= 
i 


FPZsS El == 


oe) 
I 


24] ,98 


A 
I 
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Thus 


BCE) oo= La, x(0) + ax (0) + a4x(0) + a,¥-(0)Je™ 


ie 


[b, x(0) a3 box (0) + bx(0) fy byv-(0)Je""* 


a 


Lc, x(0) + C5x(0) + C.4x(0) + cq¥-(0) JeP* (cos qt - Pa Sin qt) 


+ [d)x(0) + dpx(0) + dgx(0) + dyv,(0)JeP*sin qt (3.22) 
“¢ & = 12a lere 
with a, = 9.06x10 1 

ay = 3, 9,107" by = ~4.96x107° 

a, = 1. cea’ b, = 2.49x1077 
4 ; -7 

ay = -2,06x10 by 250 

c) = 9.78x107! d, = 1.18 

C= <2) eyquane dy = 3.93x107° 

Ca = ee dq = 8 Oe 

Cy 2.06x107" dy = 236x107" 


2. Optimization of V-(0) 


The general expression for the Integral Squared Error is 
1 


y (R-c)° dt 22) 


ISE= J 
tO 


that in this particular case, with R=0 and C =x , reduces to 


it 
ise = s xe(t) at 
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To obtain a relative minimum of this function with respect to 
V0 (0), its partial derivative with respect to the variable is made equal 


to zero, therefore 


t 
ISE : y aye) _ 
AO) a 2 x(t) nv (0 dt = 0 (3.24) 


Thus, using equations (3.22) and (3.23) 


gi 
BUOY 7! x(t)(age"" + bse" + cgeP (cos qt - pqsin qt) 


sole. 
+ dye PY sin apd: =O 
If te is taken large enough, say of the order of the settling 


time, the equation for V- (0) can be greatly simplified, since the expo- 


nential terms tend to zero. Hence the integration gives: 


ve(0) = 25865.3(a)x(0) + a5x(0) + 2,X(0)) 
-14363.1(b,x(0) + byx(0) + b3i(0)) 
-29476 .3(c,x(0) + c5x(0) + c,x(0)) 


Taking typical values of initial conditions as those obtained 


from the computer simulation, viz.: x(0) = ~6.35x107° 
x(0) = .27 
X(0) = -190 


the optimum value of v- (0) is 79.3 volts, which corresponds to an 
initial condition of -79.3 volts. It will be shown later that this 


result agrees with the experimental values. obtained from the computer. 
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To show that the inflection point of ISE vs v¢(0) 1S) sihienaee 


a minimum, the second derivative of the function is determined 


2 
NES Ss ED emt a alll i eP*(cos qt - P sin qt) 
ae ee 4 4 4 g 
Vo (0) 
+ dgeP*sin qt) dt 


‘ee fini 


A stmilar analysis could have been made to determine the optimum 
value of any other of the initial conditions or a combination of them. 
The next section will be dedicated to analyze how the wnole 


system is influenced by using different values of capacitor voltage. 
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IV. COMPUTER RESULTS 


A. DETERMINATION OF INITIAL CONDITIONS 

In order to have realistic valued initial conditions of the system 
when it enters the linear operational mode, a computer simulation of 
the bang-bang operation phase was accomplished. A convenient switching 
law was determined, such that for large enough errors, viz.: ten or 
more tracks, the motor was taken to around zero error with the appropri- 
ated final values of velocity and acceleration that would not produce a 
positive or negative overshoot when the loop was closed. 

It is important to note here that if the arm override enters a 
neighboring track, a new control law is encountered, and thus the system 
is driven completely out and away from its destination track. Therefore, 
it is imperative that once the head has entered the desired track zone 
7+ does not leave it, and if it does, it is only by a very small and 
insignificant distance (one mil or so), so that the inertia of the 
system would take care of the override, not allowing the head to be 
driven off. 

With this in mind, the results obtained from the discontinuous 
operation mode were tried as initial conditions on the closed loop mode. 
The width of the track considered was of 5 mils or 12.7x107° mm. 

Once a set of velocity and acceleration initial conditions were 
proved satisfactory to keep the head within the destination track, a 
rounded variation about these values was tried out in order to determine ~ 
the permissible amount of error in these parameters in entering the 


linear zone. It should be remembered that there is no accurate means 
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to determine velocity or acceleration, hence the need of having a 
possible wide range of these values at the entering point, and a 
control law that is based on position alone. 

A computer program for the discontinuous phase is presented in 
the appendix, the typical phase plane plots for this mode being shown 


onprig. 4.1. 


B. TYPICAL PHASE PLANE RESULTS OF CONTINUOUS SYSTEM 

The computer program used in this part of the experiment is further 
explained in the appendix. More details on the programming method can 
be found in Ref. [5} 

The simulation proved satisfactory to show the way the response of 
the system varies with the initial condition on the filter capacitor. 
By varying this initial condition over a wide range of values, the 
integral squared error as a function of the capacitor voltage was 
plotted, and the resulting curves are shown on Fig. 4.2. It can be 
seen from this figure that the integral squared error does minimize 
within the range of practicable voltages for the capacitor. Also, from 
the phase plane plots shown on Fig. 4.3, and the transient response on 
Fig. 4.4, the effect of the initial condition proves to be of great 
importance in obtaining the solution to the problem at hand. A zero 
initial condition turned out to be far from ideal in this case, and the 
voltage derived in the analytic discussion showed to be the best 
seluction in practice. 

To have a better idea of the transient response of the system, a 
computer run for a unit step input and zero initial conditions was made, 


the position vs. time curve being shown on Fig. 4.5. By a comparison 
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Phase plane plots for discontinuous operating mode of system. 
The three curves show different values of positioning orders; 
the largest used was 400 tracks, the smallest 20 tracks. 

Switching law used is a linear one, namely £ = INITIAL ERROR/7. 
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Fig. 4.2 Plot of integral squared error vs. filter initial condition 
for different pairs of velocity and acceleration values. 
Note that the optimum value for v0 (0) with the given x(t) 
and x(t) is around -60 volts. 








pig. 4.3 E vs E for x(0)=-160 m/sec and x(O)=.12 m/sec. 
Curve (a) shows response for V-(0)=-90 volts, 


(b) for V0 (0)=-85 volts, and 
((e) Rolle V- (0)=-80 volts. 
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Fig. 4.4 E vs t plot for curyes shown on Fig, 4.3. Curves (a) 
(b) and (c) show the transient response for vc(0)=-90, 
-85, and -80 volts respectively. Curve (d) has been 
added to show the transient response with zero velocity 
and acceleration, vc(0)=0 in this case. 
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Fig. 4.5 Step Input Response of the System. 


x scale 
y scale 


.O1 sec per inch 
.2 m per inch 


So 





~ 


of this curve with the ones obtained as solutions to the discontinuous 
problem, it can be seen that the initial conditions on velocity and 
acceleration greatly affect the behaviour of the system, making its 
stability quite critical depending on their values. To control this 
situation, the initial condition on the filter capacitor has to be set 
to an appropriate value which is dependent on the velocity and accelera- 
tion of the motor at the point of entry to the linear mode. Therefore, 
the best possible solution would be one that used a different value of 
capacitor voltage for each set of initial conditions at this point. 
Again, this cannot be done because of the lack of information on these 
values, thus leading to the use of a best fit voltage that will apply 
to a whole set of parameters. 

It has been the intention to show here, with the help of the phase 
plane plots, of how the error varies with the initial condition on the 
capacitor. The various computer runs showed that, although this voltage 
could prevent a negative overshoot of the error at the beginning of the 
trajectory, it produced a positive one later on if not chosen properly, 
thus providing an undesired solution. In this respect, the integral 
squared error was very efficient in detecting these effects, and 
although with the specified initial conditions of acceleration and 
velocity in some cases the system did overshoot and left the bounded 
region, using the initial condition on the filter as indicated by the 
minimized integral squared error, the reader head was maintained within 
the specified limits. Fig. 4.6 shows the difference of the trajectories 
and the overshoot caused by the use of a wrong voltage, and how it is 
remedied with an adequate oné, The values of the error when overriding 


were of the order of 0.7 to 0.9 mm., and were conveniently reduced to 
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Fig, 4.6 E vs E curves for x(0)=-180 m/sec® and x(0)=.10 m/sec. 
Curve (a) is for a vc(0)=-100 volts, (b) for 
vc(0)=-85 volts. Notice in this curve that although 
there is a tendency to override. because of a high speed 
of entry, the system is kept within the bounded region 
by the adequate use of vc(0)=0. Curve (c) shows an 
overshoot caused by the use of vc(0)=0. 


x scale = 5x107> m per inch y scale = .02 m/sec per inch 
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0.05 or 0.04. The integral squared error was reduced by a factor 
which ranges between 20 and 300, which means that the performance of 
the whole system is noticeably improved. 

Finally, it is important to point out from the experimental data, 
that when reducing the integral squared error, not only does the 
trajectory stay inside the bounded region, but also the settling time 


is minimum for this case. 
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V. CONCLUSION 


The special case presented here has served to show the significance 
of the initial condition of a filter when connected to a system to 
compensate it. In practice, other types of filters can " used, and 
a different type of feed encountered, but the analysis will be similar, 
and the effects of the initial condition should be carefully studied 
in order to obtain an optimized control. 

The bounded characteristics of this system are evidently no 
restrictions for application of the criterion here presented. It is 
generally desired to achieve minimum time of response, thus making 
the assumptions true for all cases. More sophisticated types of control 
can be used based on the results obtained, such as adjusting the 
initial condition according with the values of other parameters of the 
system if these are well known, or other convenient control that, without 
affecting cost or other factors, will provide a significant improvement 
on system performance. 

Thus, the main idea introduced by this work is that when using a 
compensator in a discontinuous manner, initial conditions of the system 
play an important role in its response characteristics, and should 
therefore be considered carefully to obtain the best possible solution. 
The specific control law used must be determined for each particular 


ease. 
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APPENDIX A 


Explanation of Simulation 


S/360 CSMP is a program for the simulation of continuous systems, 
hence limiting the possibility of simulating the Parmer discontinuous 
problem in one single program when a filter with initial condition is 
added to the diagram after time =0. This leads to the use of two 
separate programs in which, the initial conditions of the second are 
the final values of the first. 

In order to account for the initial conditions as specified by the 
problem at hand, some of the operations had to be made longer, as is 
the case of the three operations involved in solving the filter 
equations instead of just using the LEDLAG statement provided by CSM 
and the two separate integrations for the motor equations in place of 
the CMPXPL statement. 

The program provides an application-oriented input language that 
accepts problems expressed in the form of either an analog block 
diagram or a system of ordinary differential equations. The NOSORT 


statement in the discontinuous mode program allows the use of logical 


statements as the ones known in Fortran. 


40 


‘ 





LABEL MULTIMODE AUTOMATIC CONTROL SYSTEM 


DISCONTINUOUS OPERATION PHASE 


PARAM KM = 3432.3, ZETA = .9985, OMEGA = 26 
PARAM K = 55000, OMSQ = 72740.95, TAU = .00 
PARAM EINIC = (.02544,.00254) 
DYNAMIC 
NOS ORT 
INPUT = 24 
IF ( E.LE-EINIC/7 ) INPUT = -24 
SOLUTION TO INSTRUMENT SERVO EQUATIONS 
C2p0T = AUX1 + AUX2 + AUX3 
Sone; = UMIVeRL © on@s vane 
CDOT = INTGRL(O.,C2D0T) 
GC = INTGRL(0.,CDOT) 
AUX1 = -2 *% ZETA * OMEGA * C2D0T 
AUX2 = -OMSQ * CDOT 
AUX3 = KM * INPUT 
Fog —eDet 
—E = EINIC -C 
FINISH STATEMENT TERMINATES THE INTEGRATION FOR 
SPECIFIED VALUE OF CDOT 
FINISH CDOT = -.001 


PRINTOUT GF THE INDICATED VARIABLES 
PRINT C3D0T,C2D0T,CDOT»,C;E 
DATA FOR DRAW SUBROUTINE 
PREPARE E,EDOT 
TIMER FINTIM = 2.1, DELT = .001, PRDEL = .001 
END 
Sor 
END JOB 


4) 





LABEL MULTIMODE AUTOMATIC CONTRODSS sic 


CLOSED LOOP, LINEAR OPERATION PHASE 
PARAM KM = 3432.3, ZETA = .9985, OMEGA = 269.7 


PARAM K = 55000, OMSQ = 72740.95, TAU = .004 

INCON CO = -6.35E-05 

INCON C2D0TO = -160, CDOTO = .08, VCO=( -100,-85,;,0 ) 
DYNAMIC 


BE eee — ee 
SOLUTION TO FILTER EQUATIONS 


WO 1S THE VOLTAGE ACROSS THE CAPACITOR 
EC IS THE OUTPUT OF THE COMPENSATOR NETWORK 


VCDOT = -( 1 AU oy Cee 7, FA ee eee 
VC = INTGRL( eo 


O/T 

VCO 

ECG »= EPRIME =. Vv 
NT 


ERVOQ EQUATIONS 
2+ AUX3 
OTO,C3D0T) 

ne oe 


GMEGA * C2DO0T 


AO WNr 
— = 
Ut Te | el 
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EVALUATION OF THE INTEGRAL SQUARED ERROR 
E 


Ge Ee OE 1 Oo 
ER = INTGRL(O.,SQER) ~° 


mee nOUT OF THE INDICATED VARIABLES 
PRINT E2D0T, EDOT,E,INSQER 
FOR DRAW SUBROUTINE TO OBTAIN GRAPHICAL OUTPUT 
Vo~. COUT 
PREPARE E,EDOT 
TIMER FINTIM = 2.1, DELT = .001, PRDEL = .O0O01 
END 
Ser 
END JOB 


i 
= 
" 
D) 
x Il 


DATA 
Oe Ee 
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